Shell Model of Two-dimensional Turbulence in Polymer Solutions 
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We address the effect of polymer additives on two dimensional turbulence, an issue that was 
studied recently in experiments and direct numerical simulations. We show that the same simple 
shell model that reproduced drag reduction in three-dimensional turbulence reproduces all the re- 
ported effects in the two-dimensional case. The simplicity of the model offers a straightforward 
understanding of the all the major effects under consideration. 



The effect of soluble polymers on two-dimensional tur- 
bulent flows appear quite surprising Q, Instead of 
reducing the drag as in three-dimensional turbulence 
(3, 0] (thus increasing the large scale velocity compo- 
nents U 0, Q] ) , in two dimensions polymers appear to 
suppress the large scale velocity components. In ad- 
dition, the polymers affect the probability distribution 
functions of the velocity field, changing them from sub- 
gaussian to super-gaussian, with approximately exponen- 
tial tails. These phenomena appear generic, and were 
observed both in experiments in fast flowing soap films 
Q and in direct numerical simulation of the two dimen- 
sional viscoelastic equations |2(. In this Letter we argue 
that these effects can be understood using a simple shell 
model of viscoelastic turbulence. 

The shell model used to describe the effects of polymers 
on the turbulent velocity field had been derived [a, L| on 
the basis of the FENE-P model of visco-elastic flows Q . 
It was employed recently to understand successfully drag 
reduction in homogeneous three dimensional visco-elastic 
turbulence la, Ml. It reads 
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In these equations u n and B n stand for the Fourier am- 
plitudes u(k n ) and B(k n ) of the two respective fields, but 
as usual in shell models we take n = 0, 1, 2, . . . and the 
wavevectors are limited to the set k n = 2™. The nonlin- 
ear interaction terms take the explicit form 



$ n (u, B) = k n \(l - b)u n+2 B* n+l + (2 + b)u* n+l B n+2 
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with an obvious simplification for <fr„(u, u) and & n (B, B). 
Here b is a parameter which can be used to distinguish 



between three-dimensional and two-dimensional behav- 
ior (and see below for details). In accordance with the 
generalized energy of the FENE-P model, the non linear 
terms in our shell model conserves the total energy: 



E = 



n \ n / 



(3) 



The dissipative term ^fu n contains both a hyper viscosity 
at the small scales and a friction term at the large scales, 
to mimic wall friction in 2-dimensional turbulence. With 
Up = the first of Eqs. ^ reduces to the well-studied 
Sabra model of Newtonian turbulence 0- As in the 
FENE-P case we consider v v jv to be the concentration 
of the polymer c. For i/ p / we refer to the model as the 
SabraP shell model. The forcing in Eqs. is chosen to 
input a fixed amount of energy per unit time, i.e. 

K = S nJ F/u* n , F = 10- 3 V2 for / = 10, 11 . (4) 

The Sabra model with the parameter b in the range 
— 1 < b < agrees with the 3-dimensional Navier-Stokes 
equations in the sense that the energy cascade is nor- 
mal, from the large scales to the smaller scales. For 
b < b c = — 1 — 2~ 2 / 3 the situation changes qualita- 
tively, since the energy flux changes direction, going from 
smaller to larger scales, as in 2-dimensional turbulence 
[Icj . The role of the direct energy flux is taken by an 
enstrophy-like conserved variable H which cascades from 
larger to smaller scales: 
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(5) 



Thus, if the forcing is applied to an intermediate level 
kf, there are two scaling regimes with an inverse energy 
flux supported on scales k n < kf and a direct enstrophy- 
like flux for k > kf. The spectrum ^(fcn) = has 

a Kolmogorov exponents^ (k n ) ~ k n 2 ^ 3 in the inverse 
cascade regime, whereas in the direct cascade regime the 
spectrum reads (up to intermittency corrections) [l(i| 



S2(k n ) 



-2[l+log 2 (-l/l+fc)]/3 



(6) 
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FIG. 1: Power spectra &(fcn) of the SabraP model (times) 
and the Sabra model (triangles) for v v = 10~ 6 and r = 10 6 . 
Also shown is the spectrum of the B n field, i.e. (|-B n | 2 ), which 
is peaked at the Lumley scale k c . 
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FIG. 2: Probability distribution function of 3J(«i5) for the 
Sabra (squares) and the SabraP (triangles) models. The forc- 
ing is shown in Eq. @. 



When the parameter b is in the range — 1 < b < the 
coupling to the -B-field results in drag reduction, as had 
been discussed in refs. 0, 0] • We demonstrate here that 
for b < —1 — 2 -2 / 3 the coupling to the -B-field suppresses 
the large scale velocity components. In Fig. ^ we show 
the respective spectra of the Sabra and the SabraP mod- 
els, in which the forcing is identical. Note that we force 
at shells 10 and 11, and therefore we observe only the 
direct cascade, in agreement with the experimental and 
simulational works [lj, |2j . The inverse cascade regime is 
destroyed by the large scale friction and is not observ- 
able here, as in [IS]. Eq. ® predicts a spectral slope 
of -0.88 for the present value of b, while the observed 
slope is -0.96±0.05; the difference is attributed to the 
usual intermittency correction. The main point is that 
the spectrum with the coupling to the B field exhibits 
a strong suppression of the large scales, and the slope 
reduces to about -0.66±0.05. We also show, in the same 
plot, the spectrum of the B n field, (\B n \ 2 ). This spec- 
trum has a typical scale k c which is known as the Lumley 
scale IrH. 
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For wavevector smaller than k c the B n -spectrum is 
strongly suppressed, whereas for k n > k c the spectrum is 
a power law. 

The advantage of the present simple model is that the 
drag enhancement observed in the spectra lends itself 
to a straightforward interpretation. The standard new- 
tonian energy balance equation is now changed due to 
the second term in the generalized energy Q which con- 
tributes to the dissipation a positive definite term of the 
form (v p /T 2 )P 2 (B)J2 n \B n \ 2 . Thus the energy balance 



equation reads 
F = e+(v p /T 2 )P 2 (B)J2\B n \ 2 +vY, k n 4 \ u n\ 2 , (8) 



where 



(9) 



is the dissipation due to hyper viscosity. In 3-dimensional 
turbulence the dissipation is independent of the viscos- 
ity v (This fact is known in the jargon as "the viscous 
anomaly"). In the present case where the energy flux 
reverses direction e vanishes in the limit v — > 0, and for 
large Reynolds numbers can be neglected ^2 • Observing 
the fact that the last sum in © is strongly dominated 
by the forcing scale kf we can safely estimate 
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(10) 



where C is a constant of the order of unity. The first term 
in the parenthesis represents the Sabra value of \uf\ 2 . 
The subtraction of the positive definite second term is 
the explanation of the observed spectra in Fig. 

Another phenomenon that was discovered in the exper- 
iments and the direct numerical simulations is a signifi- 
cant change in the pdf's of the velocity fluctuations. We 
demonstrate that the same phenomenon is recaptured in 
our SabraP model, showing that it is generic to the in- 
teraction with the B field. In Fig. we present the 
probability distribution functions (pdf) of u n for n = 15, 
which is within the bulk of the direct enstrophy cascade. 
The qualitative change in the pdf is obvious. We quantify 
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FIG. 3: Probability distribution function of 5R(ttio) with 
Gaussian forcing on level 10, 11. Symbols are as in Fig. [5] 
Fitting a stretched exponential tails (cf. Eq. 11 U we find 
exponents 2 ± 0.15 and 1.59 ± 0.2 for the Sabra and SabraP 
models respectively. 

the difference by fitting the tails of the pdf 's to stretched 
exponential forms 

pMO] ~ exp(-/?[5R( U „)]- Q ") for R(u n ) > o n , (11) 

where a n = y (|u„| 2 ). We find ais ~ 3.45 ± 0.15 and 
ai5 ~ 1.85 ± 0.2 for the Sabra and SabraP models re- 
spectively. We note that in the forcing was Gaussian. 
We repeated our simulations for a Gaussian force and 
found results in correspondence with 0, see Fig. [21 The 
change in pdf 's indicates that the coupling to the £>-field 
results in increased intermittency. This can be quanti- 
fied by measuring the scaling exponents of the standard 
structure functions S q (k n ), 

s q (k n ) ee (KH ~ . (12) 

Indeed, the scaling exponents of the SabraP models are 
significantly more nonlinear (as a function of q) then the 
corresponding exponents of the Sabra model, cf, Fig. 0] 

The simplicity of our model again allows us to offer a 
qualitative understanding of the increased intermittency. 
Focus on the Lumley scale k c , Eq. [7| We denote u c = 
u(k c ) and B c = B(k c ). The total "energy" of the B field 
will be denoted by Q where Q = J2 n \-B n \ 2 - In light of 
the spectrum shown in Fig. \I\we can estimate Q « \B C \ 2 
up to coefficients of the order of unity. The flux of energy 
from the flow to the polymer is proportional to k c Qu c - cf. 
Eq. |21 Therefore we can write an approximate equation, 
in which we disregard the difference between complex 
numbers and their moduli: 

— w h Qk c u c . (13) 



FIG. 4: Scaling exponents ( q for the Sabra (solid line) and 
SabraP (dashed line) models. The increased intermittency in 
the SabraP model is exhibited by the increased nonlinearity 
of f 9 as a function of q. 
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FIG. 5: Probability distribution function of the quantity Q, 
see text for details. 

From this equation we see that, because k c u c w 1/r, 
that the fluctuations of Q are expected to be very strong. 
Moreover, even if u c in Eq. I)13|l were a Gaussian pro- 
cess, Q would turn out to be log-normally distributed, 
due to the multiplicative role of u c in Eq. I|13|) . Indeed, 
a direct measurement of the distribution of Q, see Fig. 
|5j confirms the existence of the anomalously long tail in 
this pdf. Finally, due to Eq. I|10(l , we see that an anoma- 
lous tail in Q would directly induce an anomalous tail 
in the pdf of Uf, as seen in the simulations. Finally, it 
is interesting to ask whether in 2-dimensional flows we 
can have enstrophy drag reduction, since this quadratic 
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invariant replaces the energy as the quantity cascaded to 
smaller scales. In our model this is impossible; our forc- 
ing with constant energy flux also produces a constant 
enstrophy flux, and as can be seen from Eq. J5J, the 
enstrophy is dominated by the larger scales (small k n ). 
Thus the reduction in energy is accompanied with the 
reduction of enstrophy. To switch the role of the scales 
to get the enstrophy to be dominated by the small scales 
(large k n ) requires changing b to compensate for the de- 
crease of \u n \ 2 with increasing k n . It is easy to check 
however that the change of roles occurs precisely at the 
value of b where the inverse cascade of energy disappears, 
i.e. b — — 1 — 2~ 2 / 3 . We thus cannot have enstrophy drag 
reduction in this model, and in our opinion neither in the 
2-dimensional Navier-Stokes case. 

In summary, we conclude that the results observed for 
the effects of dilute polymers on 2-dimensional turbu- 
lence, as observed in experiments and direct numerical 
simulations, can be readily understood with a simple 
model of viscoelastic flows once the cascade of energy is 
inverted. The same model exhibits drag reduction when 
the parameters allow a direct energy cascade. These ob- 
servations lend further support for the use of shell models 
as simple and transparent tools for understanding a va- 
riety of turbulent phenomena that remain more obscure 
when the full hydrodynamics equations are employed. 
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